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We study the statistical nature of primordial fluctuations from an anisotropic inflation which is 
realized by a vector fleld coupled to an inflaton. We flnd a suitable gauge, which we call the canon- 
ical gauge, for anisotropic inflation by generalizing the flat slicing gauge in conventional isotropic 
inflation. Using the canonical gauge, we reveal the structure of the couplings between curvature 
perturbations, vector waves, and gravitational waves. We identify two sources of anisotropy, i.e. the 
anisotropy due to the anisotropic expansion of the universe and that due to the anisotropic cou- 
, plings among variables. It turns out that the latter effect is dominant. Since the coupling between 

the curvature perturbations and vector waves is the strongest one, the statistical anisotropy in the 
5-H ' curvature perturbations is larger than that in gravitational waves. We find the cross correlation 

between the curvature perturbations and gravitational waves which never occurs in conventional in- 
flation. We also find the linear polarization of gravitational waves. Finally, we discuss cosmological 
implication of our results. 
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O . I. INTRODUCTION 

Oh' The primordial fluctuations from inflation is supposed to be statistically isotropic, Gaussian, and scale invariant. 
The nature of fluctuations is associated with the nature of de Sitter spacetime. However, since the expansion during 
inflation is not exactly de Sitter, the power spectrum is slightly tilted by the order of the slow roll parameter [l[ which 
^ ' characterizes the deviation of the expansion from the exact de Sitter expansion. The deviation from the Gaussianity is 
Ci , also known to be related to the slow roll parameter 0- On the other hand, the statistical isotropy has been regarded 
as a robust prediction so far because the cosmic no- hair conjecture is thought to be robust 0]. 
, From an observational point of view, there are various indications that there exists statistical anisotropy in the cos- 
^ ' mic microwave background radiation (CMB) Although the statistical significance of these anomalies is still under 
\ debate, the possibility of the statistical anisotropy certainly deserves further theoretical investigation Q. Recently, 
•/^ • breaking the statistical isotropy through the vector fields in an inflationary universe is proposed in the paper and 
I extended in various ways 043- However, if the vector field is relevant to inflation, it may also produce anisotropy in 
an inflationary universe whatever small it is, which seems to contradict the cosmic no-hair conjecture. 

■ From the above perspective, it is interesting to ask if it is possible to have anisotropic inflationary universe jlO MlTf . If 
, possible, it provides a simple mechanism to break the statistical isotropy by breaking the isotropy of the spacetimefisj . 

^-p ■ In the light of no-hair conjecture Q, one may deny this possibility. In fact, many attempts to construct anisotropic 
inflationary models suffer from the instability [l9| . However, recently, stable anisotropic inflationary models are found 
for the flrst time [13, [U. This can be regarded as a counter example to the cosmic no- hair conjecture. The interesting 
point is that the deviation from isotropy is related to the slow roll parameter, namely, the deviation from the exact 
deSitter expansion. Of course, that means the degree of the anisotropy is quite small. From the point of view of 

■ precision cosmology, however, it is worth exploring theoretical flne structure in an inflationary scenario. 
In this paper, we study cosmological perturbations in an anisotropic inflationary scenario we have found. The 

expected phenomenology of the anisotropic inflation is as follows: 

• There should be statistical anisotropy in curvature perturbations. 

• There should be statistical anisotropy in gravitational waves. 

• There should exist the cross correlation between curvature perturbations and gravitational waves. 

• There should be linear polarization of gravitational waves. 
The first item will be tested by the PLANCK [22|. The second one may be detected through B-mode polarization 



in the CMB [23|. The third one will imply T-B correlation in CMB l2J]. The last one could be important for the 
future direct measurement of gravitational waves through the interferometer [25j . The purpose of this paper is to 
calculate the above quantities numerically and analytically and reveal the physics behind them. Since the spacetime 
is anisotropic, the formalism treating perturbations is non-standard. Although there are many works treating the 
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cosmological perturbations in an anisotropic universe [26l - l31| . there have been several obstructions in extracting 
concrete predictions for CMB. The main obstruction was the lack of the concrete anisotropic cosmological models. 
Now, since we have such models, we have succeeded in obtaining concrete results by utilizing the canonical gauge 
which is a generalization of the flat slicing in the conventional isotropic inflationary scenario. 

Recently, during our slow preparation of this paper, two papers have appeared on the archive (3^.[33j. The first one 
[3^ studied the primordial perturbations in an anisotropic inflationary universe using a perturbative method. The 
second one [3^ investigated the same issue numerically. The conclusion is quite similar to ours. The main difference 
is the gauge used in analysis. Our canonical gauge allows us to reveal the nature of primordial fiuctuations from 
anisotropic inflation in a transparent way. Interestingly, on the contrary to a naive expectation, all of these works 
including ours imply that even if the anisotropy of the universe is very small, a large statistical anisotropy in the 
spectrum of curvature perturbations could be created. 

The organization of this paper is as follows: In section II, we review an anisotropic inflation which is caused by 
the inflaton coupled to the vector field. Here, we will see the anisotropy is determined by the slow roll parameter. 
In section III, we choose the canonical gauge which is a generalization of the fiat slicing in the conventional inflation 
and classify perturbations in anisotropic universe based on the 2-dimensional rotation symmetry. Then, we obtain 
the quadratic action for perturbed quantities. In section IV, we reduce the action to that for physical variables 
from which we can read off the structure of couplings between those variables. Based on the reduced action, we 
calculate various statistical quantities numerically and analytically to reveal the nature of primordial fluctuations in 
anisotropic inflation. In section V, we discuss cosmological implication of our results. The final section is devoted 
to the conclusion. In the Appendix A, we provide a detailed derivation of the action for 2-dimcnsional scalar sector 
perturbations. 



II. REVIEW OF ANISOTROPIC INFLATION 

In this section, we review background solutions proposed in [l^l, and see how the anisotropic inflation is realized. 
We consider the vector field Afj_ whose kinetic term is coupled to the infiaton field (f>. We note that this kind of 
model is quite natural in the context of the supergravity (sj]. The action is given by 



S = 



-9 



1 



I 



(1) 



where is the reduced gravitational constant, g is the determinant of the metric, R is the Ricci scalar, V{4)) is 
the infiaton potential, F^y is the field strength of the vector field defined by F^„ = d^A^ — d^Af^, and /((^) is a 
coupling function of the vector field. We assume that the background spacetime is given by the Bianchi type I metric 
ds^ = —dt'^ + ^ii^)^ dx^^ . As for the vector field, it can be reduced into "electric" (Foi) and " magnetic" (F^ ) 

components, and here we consider only the "electric" component for simplicity. It is not hard to prove that the 
direction of the "electric" field docs not change in time by solving its evolution equation. Without loosing the 
generality, one can take x-axis in the direction of the "electric" field. Using the gauge invariance, we can express the 
vector field as A^j^dx^ = v(t)dx. Thus, there exists the rotational symmetry in the y-z plane. Given this configuration, 
it is convenient to parameterize the metric as follows: 



(2) 



where e", a and (t_ are an isotropic scale factor and spatial shears, respectively. We also define the averaged Hubble 
parameter as H = a. Here, the lapse function J\f is introduced to obtain the Hamiltonian constraint. With the above 
ansatz, the action becomes 



S 



(3) 



where an overdot denotes the derivative with respect to the physical time t. First, its variation with respect to (t_ 
yields 



(4) 



This gives &- (x e hence, the anisotropy in the y-z plane rapidly decays as the universe expands. Hereafter, for 
simplicity, we assume cr_ = and set the metric to be 



ds^ 



(5) 
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Next, the equation of motion for v is easily solved as 



-2„-Q-4o-„ 



PA , 



(6) 



where pA is a constant. Taking the variation of the action with respect to Af, a, a and (f) and substituting the solution 
([6|), we obtain the following basic equations: 
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\-2„-4q-4<t 



-3da + ^/(0)-2e-'^"-4- 



(7) 

(8) 

(9) 
(10) 



where the subscript in denotes derivative with respect to (f). Let us check whether inflation occurs in this model. 
Using Eqs. (O and (jS]), the equation for acceleration of the cosmic expansion is given by 



(e«)- 
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(11) 



We see that the potential energy of the inflaton needs to be dominant and the energy density of the vector field 
Py = p^/((/))^^e^*"^*°'/2 and the shear S = tr should be subdominant for inflation to occur. Next we want to see 
if the anisotropy is produced during inflation. Here we look at the ratio of the shear to the expansion rate "E/H to 
characterize the anisotropy of the inflationary universe. Then, Eq.([9]) reads 



-SHE 



2^2 
~3 
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If the anisotropy converges to a value, i.e. S becomes negligible, the terminal value should be given by 

S 2, 



H 



:7^ 



(12) 



(13) 



where we used the slow roll equation ~ K'V{(t))/3 derived from Eq.(l7]) and defined the energy ratio TZ = pv/V{(j)). 

In order to realize the above situation, py must be almost constant. Assuming that the vector field is subdominant 
in the evolution equation of the inflaton field Eq. pH)) and the conventional single field slow- roll inflation is realized, 
one can show the coupling function f{(f>) should be proportional to e"^" to keep p^ almost constant. In the slow roll 
phase, e-folding number a is related to the inflaton field (p as da — ~K''V{(j))d(j)/V^ as usual. Then, the functional 
form of f{<j)) is determined as 



2.' J ^0 



(14) 



For the polynomial potential V oc (/>", for example, we have / = e" 



from the condition that the vector field is negligible in Eq. pH]) . that is, \p\f{4>) 



In this scenario, the anisotropy is restricted 



-3 , 



< IV^I. Substituting 



Eq. ([H)) into this, we obtain TZ <^ ey/2, where ey is the slow- roll parameter defined by ey = 1/2k^ (V^/V) . Thus, 
the anisotropy is constrained by E/H <^ 0{ev)- 

The above case is, in a sense, a critical one, and next we want to consider beyond it. For simplicity, we parameterize 
fW by 

/(,/,) = 6'""'^ (15) 
where c is a constant parameter. Now, we look at what happens when c > 1. The basic equations become 



= a + a' = -2(7 
— — 3d(/) — 



3 ^ 

1 n 



(16) 
(17) 
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In this case, if the vector field is initiahy small TZ <C ey/2c, then the conventional single field slow- roll inflation 
is realized. During this stage / (x e"^"^" and the vector field grows as p„ oc e^(<^-'^)°'. Therefore, the vector field 
eventually becomes relevant to the inflaton dynamics Eg. (fT7|) . Nevertheless, the accelerating expansion of the universe 
will continue. The point is that TZ cannot exceed ey/2c. In fact, if TZ exceeds ey/2c, the inflaton field (j) does not 
roll down, which makes p„ = p\f{(j>)^'^e~'^°'~'^"'/2 decrease. Hence, ^ y{4') always holds. In this sense, there 
exists an attractor where the inflation continues even when the vector field affects the infiaton evolution. The inflaton 
dynamics is determined by solving the slow-roll equation: 

- 3a<P ~V^+ pi/- V^e-^"-^^ = 0. (18) 

Assuming = K^V{(f>)/3, this yields 

d0 ^ ^ _ ^ 2c^e''"~'"~''"' ^ ^''^ (19) 
da a k'^V 

This can be integrated by neglecting the evolutions of 1/, V^, cr as 

l + const.e-^(^-i)"+^" 



' " c-1 

Substituting this into the slow-roll equation Eq. ([T9| , we obtain 

d(f> _ V^c-lV^ 



da kW c kW 



1 + const. e" 



-4(c— l)ct+4cr 



(20) 



(21) 



This clearly shows a transition from the conventional single field slow-roll inflationary phase, where dcp/da = —V^/k^V 
holds, to what we refer to as the second inflationary phase, where the vector field is relevant to the inflaton dynamics 
and the inflaton gets 1/c times slower as d(/)/da = —Vtp/cK^V. In the second inflationary phase, the energy density of 
the vector flcld becomes 

P. = fe— = (22) 

which yields the anisotropy of E/iJ = 27^/3 = (c — l)ev/3c^. From Eqs.Q and ([8]), the slow-roll parameter defined 
in terms of the scale factor becomes 

- " 1 1 2 ;2 2 2 \ 1 
£// = — = "o'* - P^l = -^v , (23 

a^ a^ \ 2 3 / c 

where we neglected the anisotropy and used relations 4)/a = d(t)/da = —V^/cn^V and dP' = k^F((/))/3. Thus we have 
a remarkable result 

In the next section, we will make a perturbativc analysis during this second inflationary phase. 



III. CANONICAL GAUGE IN ANISOTROPIC INFLATION 



Since the background is expanding anisotropically, we cannot use the standard cosmological perturbation theory. In 
this section, we classify perturbations under the 2-dimensional rotational symmetry and obtain the quadratic actions 
for 2-dimcnsional scalar and vector sectors. In order to grasp the meaning of variables, we start with the isotropic case 
and make a gauge transformation from the flat slicing gauge to the appropriate gauge for 2-dimensional classification. 
Then, the resultant gauge can be promoted to the anisotropic spacetime. The gauge we have chosen makes the 
analysis and the interpretation easier. Once the gauge is fixed, the quadratic action can be calculated. 



A. Gauge Fixing and Classification of perturbations 



First, we will start with the spatially homogeneous and isotropic universe. For simplicity, we consider flat space. 

ds^ = a^{ri) [-dif + S^jdx'dx^] , (25) 



5 



where we took a conformal time 77. In that case, we can use S-dimensional rotational symmetry to classify the 
perturbed metric. When we want to have diagonal quadratic action, we take the following gauge 



ds^ = [-(1 + 2A)d'n^ + 2{B^i + V^)dridx^ + (Sij + hij)dx'dx^] 



(26) 



where wc imposed Vi,i = and hijj = ha ~ 0. If we ignore vector and tensor perturbations Vi, hij, the above gauge 
is called the flat slicing gauge. Now, let us move on to the Fourier space. Since there exists 3-dimensional rotation 
symmetry, we can take a wave number vector to be k = (A:,0,0). Then, the perturbed metric has the following 
components: 



/ -2a^A 




a^V2 















* 








V * 









* is symmetric part. 



(27) 



Here, we utilized the special choice k = (fc,0, 0) to solve the constraints Vi^i = and h^j j = = 0. With the 
same reason, only remains. We defined hyz = —hx, hyy ~ —h^z = Now, we will pretend that we have only 
2-dimensional rotation symmetry in y-z plane. In that case, at best, we can take k = (kx,ky,0). Hence, we make a 
rotation in the x — y-plane so that the wave number becomes k = (fc^,, ky,0). 



kx 
ki, 



kx 





(28) 



where we have a relation fc^ = k1 



ky. Under this rotation, the perturbed metric becomes 



/ -20^ A f^a'B^-^a^Vo 



V 



O ■p-'l_|_ 



2 «T I, 



a^!fhx 



(29) 



To simplify the perturbations, we can make use of gauge transformation 



(30) 



where the semicolon denotes the covariant derivative with respect to the background metric. Taking the parameter 



kx 
2iP 



2k^' 



we obtain 



/ -2a^A ^a^B,x-!fa^V, 



\ 



a h+ 
* 
* 



2ik 



k " 



0?h4 



2ik? 







J 



(31) 



It should be noted that wc did not change slicing but performed only the spatial coordinate transformation. Therefore, 
we arc still working in the flat slicing where the 3-dimcnsional scalar curvature vanishes. 

In our anisotropic inflation models, the available symmetry is actually small. The background metric is given by 



dsl = a(77)2(-d7?2 + dx^) + hi-qfidy^ + dz^), 



(32) 



that is, a 



, h = e""'"'^, drj = dt/a. Notice that the conformal time in anisotropic inflation is the conformal time 



in 2-dimensional part {r],x). Even in this anisotropic spacetime, as we have done in pip , one can take the following 
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gauge 

(^.goo 5gQ^ Sgoy Sgo^ \ 

* Sgxx Sgxz 
* Sgyy ' ^^^> 

* * * 5gzz J 

where we can impose further conditions so that the perturbed metric goes back to pip in the isotropic hmit. 

One can classify the perturbed metric using the rotational symmetry my — z-plane. In 2-dimensional flat space, 
an arbitrary vector to" where a = y, z can be decomposed into the scalar part to"^ ^ and the vector part to"^ = 0. 
Since there exists no tensor part in 2-dimensions, 2-dimensional tensor can be constructed from the 2-dimensional 
vector. Because of the symmetry, the scalar and vector parts are not mixed in the equations. Thus, the metric 
perturbations can be classified to the scalar sector and the vector sector. Thanks to the symmetry in the y — z plane, 
without loss of generality, we can take the wave number vector to be k = {kx,ky,0). Hence, the vector sector in 
2-dimensional classification can be identified as Sgoz, Sgxz in the above perturbed metric. The remaining components 
Sgoo, Sgox, Sgoy, Sgxx, Sgyy, Sgzz belong to the scalar sector. 



Sgiiu = 

V 



B. 2d vector sector 



Thus, the perturbations that belong to 2d vector perturbations, can be written as 





V 



* 



IPV 

* / 



(34) 



where we have incorporated the anisotropy while keeping the spatial scalar curvature to be zero. As to the vector 
field, we can take 



^^vcctor 



(0 ,0 ,0 ,D) 



(35) 



Note that we have no residual gauge transformation and, in particular, D is a gauge invariant under abelian gauge 
transformations. And, as we have seen in (j3ip . T corresponds to the cross- mode polarization of gravitational waves 
in the isotropic limit a ~ b. 

Using this gauge, we can calculate the quadratic action as 



= / drjd^x 



-'3,x 



4a2 



2a2 

2&2 2-' '"^ 2 



■D'T 



(36) 



Since the perturbed shift function jS^ does not have a time derivative, it is not dynamical. There are two physical 
degrees of freedom T and D in this 2-dimensional vector sector. 



2d scalar sector 



For the 2-dimensional scalar sector, we define the metric perturbations 

/ -2a2$ aPi a/32 
* 2a^G 



V 



262G 





-2fo2G / 



(37) 



where we have kept the spatial scalar curvature vanishing. The scalar perturbation will be represented by S4>. The 
variable G and S(f> are the gauge invariant variables that correspond to the plus mode of gravitational waves and the 
scalar perturbations, respectively, in the isotropized limit a = b. And, we set the perturbed vector to be 



SA 



scalar 



(Mo ,0 , J ,0) , 



(38) 
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where we have fixed the abchan gauge by putting the longitudinal component to be zero. From these ansatz, we can 
calculate the quadratic action as 



^scalar 



d xdrj 



2a2 



.62 



1 1 bb' 6^ 1 ( a! y 



1. 



1 



,.'2jCJ,2 „2i,2t^^2 



2o? 

fv''^G' - 2a^6^F$G - 2bb'<^'G - ( ^fU^'^ + a-'b'V^, ]S(j){G + <P)+ b^(t>'5(t>' (G - $) 



6^ //,2 

-fv'^G^ - 2a^b^V^G - 2bb'^'G - 

2a^ \ 



(39) 



Here, S^'^'^^^^ consists of $, /32, G, JAq, (5(/) and J. Among them, $,/3i,/32 and ^Aq are non-dynamical and can be 
eliminated. 

IV. THE NATURE OF PRIMORDIAL FLUCTUATIONS 

Now, we are in a position to calculate the statistical properties of primordial fluctuations from anisotropic inflation. 
For this purpose, we need to reduce the action to the one for physical variables. Then, we can quantize the system 
and specify the vacuum state. We analyze the vector sector and the scalar sector separately. 

A. Action in slow roll approximation 



First, let us eliminate non-dynamical variables from the action for the 2-dimensional vector sector. In Fourier space, 
the action for 2-dimensional vector sector becomes 



/■2 , f2r,2 ^2 /■2,//,2 

.__|r'|2 _ ^fc^irp + i-\D ? - 4^fc2|^|2 _ 4fc2|^|2 _ 



4a2' 



52. 



262 



D T* 



(40) 



where we omitted the Fourier indices k for simplicity. After making the above action real manifestly and completing 
the square of /^a, we obtain 



4a2 



/fc2 



52 B 



^2 i-2 .2 



f2 , f2 f4,/2 1,2 



2 fc2 



{v' D* + V*' 



(41) 



where k is given by ^(77) — +a'^k'^/b'^ which becomes constant in the isotropic limit a — b. The first squared 
term vanishes after substituting the equation of motion for /^a. Now, we define canonically normalized variables as 



^ ^r, D ^ fD. 

V2k 



(42) 



Then, using the canonical variables and integrating by parts, we obtain the reduced action for physical variables 

(k/b) 



V2 a 5M'''*^'*''^TW^'^*^'*^) 



(43) 
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In the isotropic limit a = b, T and D represent the cross-mode of gravitational waves and vector waves, respectively. 
The second line in the action describes how both waves are interacting to each other. 

Next, we use the slow roll approximation to simplify the action. To obtain the homogeneous background metric, 
we integrate the following equations 



H 



H = 3''" 



by assuming e^/ch ^ a /a. The resultant expressions are 



(44) 



(45) 



In this approximation, the universe shows anisotropic power law inflation. We should recall, in the second inflationary 
phase, the variable / is given by 



c- 1 



(46) 



Note that the range (1,cxd) for c corresponds to (0,1) for /. The above approximation gives useful formula for the 
subsequent calculations 



a 

52 

k' 



a 

(-?7)-2 [ 1 + 2e^ + 



= (-r,)-2 [ 2 + 3eff ] , = (-r;)-^ [ 1 + 26^ + I^h ] 

a 

, ^ = (-7?)-2 [ 2 + 3eff + ] , 

h" fc2 n2 



From the background equations Eqs.((7|)- (|T0)) . it is easy to obtain 



2„,'2 



^ 6" ^ a'2 a" I) 



IP 



ah 



Using the formula (jTZ]). we obtain 



2,, 2 



3{-r])-^IeH 



(47) 



(48) 



(49) 



From Eq. ([6|), the background equation for the vector can be found as 







From this equation, it is easy to deduce the relation 



/ 



= i-r,)-^ [-2 - 3eH + m- ^I^h] , 



where r]H is a slow-roll parameter defined by 

^ = 2^ i2eH - m) = 2(2eH - 7?ff)(-7?)-i 

Of course, rjn is not related to the conformal time 1^. Furthermore, we obtain 

f" 

— = (-7?)"'[2 + 9eH-3r/ff +6/eff] . 



(50) 



(51) 



(52) 



(53) 
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Substituting these results into the action, we obtain the action in the slow roll approximation 

^vector ^ I ^^^^3^ 



^\f\^ + ^[~k^ + i-7])-^2 + -3eH + ^IeH + SIeHsmH}] jfp 



1, 



\D 



' |2 



1 



[-k^ + {-r^y^ {2 + 9eH - 3r]H + Glen sin^ 9}] \D\ 



-{~r])-'sm9{T D* +T* D) 



{~jjy^sme{TD* + r*D) 



(54) 



where we have defined sin0 = kya/kb. This 9 represents the direction dependence. In the isotropic limit 1 = 0, the 
Lagrangian for T becomes the familiar one for gravitational waves in a Friedman-Lemaitre universe. 

In a similar way, we can derive the quadratic action for physical variables in the 2-dimensional scalar sector. 
The details can be found in the Appendix A. Moreover, it is straightforward to deduce the action in the slow roll 
approximation. The resultant action is given by 



^scalar 
^GG 



\G 



' |2 



[-fc2 + |2 + ien + ilen + S/e^ sin^ 9]] |G|^ 

+ ^ [-^' + {2 + 9eff-3r;^^-6/eHsin2 0}] jjp, 



37]H 12/ 



l-I l-I 



+ 12/e 



H 



24/ 
l-I 



sm 



= -3/ 



^^{-f])-^ sin^ 9 {GS<I>* +G*6(f>) , 



^(-.)-sm.(.r./ + ^0V-)-J^(-.)-sin.(.^V-+.^J*) 



-JG 



(-7?)-^ sin 61 G* J + G J 



(-77)-2 sine (G* J + GJ*) 



where we defined canonical variables 



G = V26G , J : 



■J 



(55) 
(56) 
(57) 
(58) 

(59) 

(60) 
(61) 

(62) 



Here, G , J and S(f> represent the gravitational waves, the vector waves, and the scalar perturbations, respectively. 
The above action shows there exist the interaction among these variables. We notice the scalar part contains / 
without suppression by a slow-roll parameter en- Therefore, to obtain the quasi-scale invariant spectrum of curvature 
perturbation, / itself has to be small. 

From the actions ([54)) and ((55|) . we see there are two sources of statistical anisotropy of fluctuations. First, the 
statistical anisotropy of fluctuations comes from the anisotropic expansion itself. Intuitively, this can be understood 
from the anisotropic effective Hawking temperature Hcs/2-k, where denotes the effective expansion rate. Indeed, 
the expansion rate in the direction of the background vector is relatively small, hence the effective Hawking temper- 
ature is low. Then, this direction has less fluctuation power compared to the other directions. Thus, the effective 
Hawking temperature induces the anisotropy in the power spectrum of fluctuations. This effect is encoded in (|56p . 
([57)) . and ([55)) . The other source of the statistical anisotropy of fluctuations comes from the couplings (|5^ . ()60p and 
()6ip due to the background vector field. The essential structure of couplings can be understood without complicated 
calculations. Take a look at the following term 



Here, we should recall the order of magnitude of background quantities 



(63) 



/ 



For example, to obtain the J — G coupling, one of F^j^y have to be replaced by the background quantity v' . Hence, 
the coefficients in the J — G coupling should be proportional to fv' which is of the order of ^/Ieh- This explains the 
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strength of the couphng in (|6T|) . Similarly, J — 5(f) coupling should be proportional to frf^v' because we have to take 
the variation with respect to (j). Hence, wc can estimate its magnitude to be a/T- This agrees to the interaction term 
(pO)) . Finally, the coupling G — 6(j) has a magnitude of the order of f^v''^ which is proportional to I^JTh. This shows 
a good agreement with the coupling (j59p . Thus, we can understand why there is a hierarchy among the couplings of 
the gravitational waves, the vector waves and the scalar field. 



B. Numerical results 



In this section, wc will calculate power spectrum of various variables. To set the initial conditions, we need to 
quantize this system by promoting canonical variables to operators which satisfy the following canonical commutation 
relations 

[Q,(7?,x),Pfc(7,,y)] =i(5,bJ(x-y) , \Qa{^^,^)Mn,y)\ = \Pa(j].^).Pb(j].y)\ = ^ (64) 



where Qa (1 < a < 5) denote the variables F, Z?, (50, G, J in this order, while Pa is their conjugate momentum 
defined by Pa = 5L/5Q'^. The point is that, with a given wave number, the actions and ([55)) reduce to those of 
independent harmonic oscillators in the subhorizon limit — fcyy 3> 1 

^ = E ^ / '^^^'^ - k^W] ■ (65) 

a=l ^ ^ 

We choose the Bunch-Davis vacuum state |0) by imposing the following initial conditions at an initial time rji 

Qaim) = \/^("'^.k + al,-k) ' Q'aiVi) = -i^iaaM ~ 4,-1.) ' (66) 
where aa,k is an annihilation operator whose commutation relations are given by 

SabS^^Hk-k'), Kk,a;,,k'] =0 . (67) 



It is easy to verify the commutation relations imposed for Pa ,Qa- The Bunch-Davis vacuum |0) is defined by 

fla.klO) = 0. 

Later, the system evolves according to the actions ([54]) and ([55]) and the variables Qaiv) are given by linear 
combination of Ub^k as 

5 

Qa,k(??) = [cafc(?])a6,k + c*t(??)a[_k ' (68) 

b=l 

where the transfer matrix Cab depends on the wave number k and the time rj. The time dependence is determined 
by solving the classical equations of motion with the initial conditions given by the coefficients in Egs. lpSl) . And the 
expectation values of operators in this vacuum state are evaluated as 

5 

{0\Qau{v)Q*bp{v)\0) =Y.^-d'=*bd6{k + p) . (69) 

d=l 

Especially, we are interested in the power spectrum of the scalar perturbations, cross mode and plus mode of gravi- 
tational waves, the cross correlation between scalar perturbations and the plus mode of gravitational waves, and the 
linear polarization of gravitational waves. The power spectrum of scalar perturbations is given by 

5 

(O|,5<^k(?7)<50p(r/)|O) =^|c3d(7?)P<5(k + p) = P,-4k),5(k + p) , (70) 

d=3 

where we took into account the fact that the 2-dimcnsional vector sector is decoupled from the scalar sector. The 
power spectrum of the cross and plus mode of gravitational waves read 



(0|fk(r,)f p(,?)|0) = J2 Mv)\'Sik + p) ^ Pr(k)^(k + p) , (71) 

5 

(0|Gk(r/)Gp(ry)|0)==^|c4<j(77)P<5(k + p)=PG(k)<5(k + p) . (72) 



d=3 
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Here, we used the decoupling of the 2-dimcnsional scalar sector and the vector sector. We can also calculate the cross 
correlation between the plus mode of gravitational waves and the scalar perturbations 



(O|50k(?7)Gp(7?)|O) =^C3d(?7)c:rf(?7)<5(k + p) = n>G(k)<5(k + p) . 



(73) 



d=3 



The linear polarization of gravitational waves can be calculated once the power spectrum of the cross and plus mode 
are calculated. 




0.001 



0.0001 



1 10 

time(e— folding number) 



50 



FIG. 1: Evolution of anisotropy of curvature perturbations. Here we depicted the anisotropy Ps4>(ii-)\g_^^2/ ^^'t'(^)\g-o ~ f 
a function of time. We set e-folding number to be zero at the time of horizon crossing of the given mode. The both axes are 
taken in log scale. 




1 10 

time(e— folding number) 



FIG. 2: Evolution of anisotropy in gravitational waves. Here we depicted the anisotropy Px {k)\g_^^2/ (^)\e-o ~ ^' ^^ere 
X = r, G, as a function of time. The axes are in log scale. As one can see, the difference between two modes is quite small. 



Now, we numerically calculate the transfer matrix Cab with parameters = 10~^, / — 10^'^ and evaluate various 
statistical quantities. First of all, we need to calculate the statistical anisotropy in the curvature perturbations. We 
depict the time evolution of the anisotropy in the curvature perturbation 



e=7r/2 



(74) 



as a function of time in FigHJ Taking look at FigU we see the anisotropy grows with the squared of the e-folding 
number. The amplitude is larger than the expected one due to the coupling proportional to a/T which is larger than 
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0.1 1 10 50 

time(e— folding number) 

FIG. 3: Evolution of cross correlation between curvature perturbations and the plus mode of gravitational waves. Here we 
depicted the value —Psr/>G0^)\g_^^2/^^'l'^^'^\g—ir/2 ^ function of the time. The axes are in log scale. 



/2 



9=7r/2 




5 10 20 30 40 60 

time(e— folding number) 

FIG. 4: Evolution of the linear polarization of gravitational waves. Here we depicted the value fr(k) |g_^^2/^'G(k) |g_^^2 
as a function of time. The axes are in log scale. 



the expected one \/IeH- Similarly, in FiglJl we show the power spectrum of the cross and plus mode of gravitational 
waves 



/2 



1 , 



(75) 



8=0 



where X = F, G. Next, we are interested in the cross correlation between the curvature perturbations and the 
gravitational waves. In Figl31 we depict the cross correlation between the curvature perturbation and the plus mode 
of gravitational waves which is normalized by the power spectrum of scalar perturbations 



e=7r/2 



(76) 



/2 



In FiglSl we see the cross correlation between the scalar perturbations and the gravitational waves grows quadratically 
in the e- folding number. The amplitude is proportional to the coupling between G and (50, i.e., l^ftn- The difference 
between the power spectrum of two polarization modes normalized by the power spectrum of the plus mode of 
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gravitational waves 



/2 



(77) 



= Tl/2 



characterizes the hnear polarization of the gravitational wave. As you can see from Figl2]and FigHl the anisotropy of 
power spectrum grows quadratically in the e-folding number, while the difference between the two polarization modes 
shows the growth of higher power in e-folding number. 



C. Analytical Estimation 



In this section, we treat the anisotropy perturbatively and estimate its magnitude by using perturbation in the 
interaction picture. In the interaction picture, the expectation value for a physical quantity 0{r]) is given by 



(m 10(77)1^71) (0 



Tcxp ( i I Hj{ri )dij 

Vi 



oiv) 



Tcxp ( —7 / Hj{ij )drj 

Vi 



(78) 



where \in) is an in vacuum in the interaction picture, T and T denote a time-ordered and an anti-time-ordcrcd product 
and Hi denotes the interaction part of Hamiltonian in this picture. This is equivalent to the following 

{in\0{ij)\in)=Y,^'' dm dm^i--- dm {0\[Hji,n),[Hiiv2), ■ ■ ■ [Hi{m),o{v)] ■ ■ ■]]\o) . (79) 



The successive approximation is made by truncating the summation at a certain order A^. In our analysis, we assume 
the noninteracting part of Hamiltonian to be that of free fields in dcSittcr spacctimc 



5 r 



^IQ:p-^(^^-2(-7y)-2)|Q„P 



(80) 



and the operators in the interaction picture are given by 

Qn.Miv) = '«(77)a„,k + u(77)*a]'^ 



1 ~ikr, ( l_ 1_ 



ki] 



(81) 
(82) 



And the rest of the Lagrangian ([M)) -(|5T |) is regarded as the interaction part L/ = L^^^ — Lq. To see the leading 
effect on the anisotropy in the scalar perturbation, which is of the order of /, we evaluate the correction due to the 
interaction given by 



d'^k 



^^{-ri)-^ sine (f/ J + S~(t>' .P ] + 



(83) 



Note that in the analogy with the slow-roll parameter in the ordinary slow-roll inflation, the term proportional to 



I sin^ 96(1)6(1)' in ([58]) can be expected to give the anisotropy ^(m |(5</)k(5(/)p| m)/(0 |(5(/)k<5(/)p| 0) ^ sin^ 9IN{k) where 
N{k) is the e-folding number from the horizon exit. Thus, the leading correction comes from the interaction through 



the term Hf'\ The leading correction is given by 



6{in\6(t)i^{i-l)6(j)p{r])\in) = i'^ I di]2 diji (o Hf\rii), Hf'' {'1-12), 6(j)]^{r])6(t)p{r]) 



(84) 



Using Eqs. (|8ip and commutation relations for the creation and annihilation operators, we obtain the anisotropy 
expressed as follows 



d{in 


\6(j)kS(l)p 


in) 




6(j)k6(j)p\ 


0) 



iv) 



24/ 

1 - / 

xlm 



dri2 / drji- 

Vi •'Vi 



-Im 



-(-772) {■q2)u*{-q) + {-ri2) '^u{r]2)u* {7]) 



l(?7l)u*(?72) |-(-7?l) {t]i)u*{7]) + {-TJi) '^u{tJi)u*{tj)^ 



(85) 
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where Im denotes the imaginary part. Substituting the function form of u (j82[) and introducing time variables x = kr], 
Xi = kiji and X2 = kr/2, we have 



5{in 


\S(t)kS(j)p 


in) 






0) 



6/ 



sin^6i 



dX2 



dxi 



1 1 



1 



cos(-2xi +X + X2) ( 1 + ^ — 



-sin(-2xi + X + X2) - 



-Xi -X2 
1 

XXi XX2 X1X2 

1 111 
\ 

XX1X2 Xi X X2 



cos(-X2 + x) - sin(-X2 + x) 



(86) 



The contribution to the integral from the subhorizon — xi 3> 1 is negligible. In the limit of superhorizon —xi <C 1, we 
also have — X2 <C 1,— X ^ 1- Hence, the integrand in Eq. (|86p approximately becomes 8/xiX2- Thus, the anisotropy 
can be evaluated as 



5 {in 


\54'\^5(f)p 


in) 




6(i>\^5(i)p\ 


0) 



(x) 



6/ .2. 
: sm 



1 -/ 
24/ 



X 



1 



C^X2 

sin^e* N'\k), 



^2/1 Ar2 / 



X2 8 

«Xi 

-1 X1X2 



(87) 



where N{k) = — ln(— fcr;) is the e-folding number from the horizon exit. 

For the anisotropy of both two polarizations of gravitational waves, the similar calculations give 



Men sm^ N'^{k) , 



where we used the interaction term in the action ([M)) for F and that in (pT|) for G. It is interesting to calculate the 
cross correlation. The leading contribution comes from iJ/^ and uf. The result is as follows: 



S{in 


fkfp 


in) 5 {in 


GkGp 


in) 


(0|fkfp| 


0) (0|GkGp| 


0) 



{in 


(50k Gp 


in) 


(0 


(J^k-^^pIO) 



^24/ 



1 



iV^(fc) 



(89) 



As we will soon see, this might give a detectable number. On the other hand, it turns out that the linear polarization 
of the gravitational wave must be the effect of the order higher than /. 



V. COSMOLOGICAL IMPLICATION 



Now, we are in a position to discuss cosmological implication of an anisotropic inflationary scenario. As we have 
listed up in section I, there are many interesting phenomenology in the anisotropic inflation. Here, we summarize our 
results: 

• There exists statistical anisotropy in curvature perturbations of the order of 2AIN'^{k). 

• There exists statistical anisotropy in gravitational waves of the order of QlenN^ik). 

• These exists the cross correlation between curvature perturbations and gravitational waves of the order of 
~2M^N-'{k). 

• There is small linear polarization of gravitational waves. 

Due to the interaction on superhorizon scales, there is an enhancement factor N'^(k) in the above quantities. Because 
of this enhancement, even when the anisotropy of the spacetime is quite small, say E/iJ ^ 10"'' in our example, the 
statistical anisotropy imprinted in primordial fluctuations can not be negligible in precision cosmology. 
Let us be more precise. The anisotropy is often parameterized by 

P(k) = P{k) [1 + sin^ 61] . (90) 

We have estimated 5* ~ 2AIN'^{k) which could be large. The current observational limit of the statistical anisotropy 
for the curvature perturbations is characterized by < 0.3 [l^. Let us suppose = 0.2. Then, we have the 
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anisotropy in the gravitational waves with g^, ~ 10~^. The hnear polarization is of the order of 10~^. Using the 
definition of curvature perturbations C = Scj)/ one can translate the cross correlation (|89| between the scalar 
perturbations and gravitational waves to that between the curvature perturbations and gravitational waves normalized 
by the power spectrum of curvature perturbations: 

\ 1'-'^ PI ' ~ _20/Ar2(fc)g^ _ ^ _2 X 10-^ , (91) 



(o CkCp 0) 



where we used ^ 0.2 and ejj ~ 10~^. Since the current constraints on the TBjTE ratio is of the order of 10"'^ [l[, 
we need to improve the accuracy by one more order, which might be achieved by PLANCK. 

In [3^, it is pointed out that the sign of is different from the observed one. However, it might be possible 
to modify the model so that the sign of is flipped. For example, we can consider two vector fields. Then, the 
orthogonal direction to the plane determined by two vectors becomes a preferred direction. In this case, we can expect 
the sign of 5* becomes opposite. We can also utilize anti-symmetric tensor fields to achieve the same aim. The details 
will be reported elsewhere. 



VI. CONCLUSION 



We have studied the statistical nature of primordial fluctuations from an anisotropic inflation which is realized 
by a vector field coupled to an inflaton. First, we have classified metric fluctuations according to the 2-dimensional 
rotational symmetry. To choose a convenient gauge in an anisotropic universe, we have started from the flat slicing 
gauge in an isotropic universe and made an appropriate gauge transformation to get a canonical gauge in an anisotropic 
universe. This gauge choice has made the subsequent analysis and the interpretation of the variables easier. Using the 
canonical gauge, we have revealed the structure of the couplings between curvature perturbations, vector waves, and 
gravitational waves. We found that there are two sources for anisotropy, i.e. the anisotropy due to the anisotropic 
expansion of the universe and that due to the anisotropic couplings among variables. It turned out that the latter 
effect is dominant. We have numerically obtained power spectra. We also presented analytical formula using in- 
in formalism. Since the coupling between the curvature perturbations and vector waves is the strongest one. the 
anisotropy in the curvature perturbations is larger than that in gravitational waves. More interestingly, we found the 
cross correlation between curvature perturbations and gravitational waves which is peculiar to anisotropic inflation. 
We also found the linear polarization of gravitational waves. Although there are several mechanism to produce circular 
polarization in the primordial gravitational waves [ssj . this is the first example which realized the linear polarization 
in the primordial gravitational waves. 

We have only considered power spectrum for simplicity. However, as is pointed out in the paper @, the statistical 
anisotropy could appear in the non-Gaussianity strongly and modify the shape of the bispcctrum and trispcctrum. 
Hence, it is interesting to study non-Gaussianity in anisotropic inflation models. 

We can extend anisotropic inflation in various ways. Although we have investigated a chaotic inflation in this paper, 
it is easy to extend the analysis to other inflation models. It is possible to incorporate multi-vector fields. From the 
string theory point of vieWjji is intriguing to consider anti-symmetric tensor field. It is also interesting to consider 
other Bianchi type models |36| in the context of anisotropic inflation. 
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Appendix A: derivation of reduced action for 2d scalar sector 



In this appendix, we present a detailed derivation of the action for the 2-dimcnsional scalar sector. Here, we need 
to eliminate non-dynamical fields (5Ao, fii and $. 
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First, let us see the terms concerning SAq in Fourier space 



2a2 



-ik 



By completing the square as 



(Al) 



2a2 



6An 



ikj-v 
~fc2~ 



/ 



.ky a I 



2a2 fc2 



fc^ru G + $ - 2 



-fc / 



(A2) 



the variable 5Aq can be eliminated because the first squared term vanishes after substituting the equation of motion 
for 6Af). Similarly, the terms containing /3i are given by 



kim^+Pi 



1 6V' 
'-k^kyl32 + ikx — $ - ikx——( 
4 a 2a 



c.c. 



(A3) 



Completing the square gives 



kl 



k'x , kxbb ^ „ kx b'^<t> 

Qi Ky Qj 



kx 
4 



^2 - $ + 2i- ^ 



(A4) 



Again, we can ignore the first term. Taking a look at the terms related to the variable (32 , we see only the linear term 
in /32 appears as 



/32 



b 9 k 



e 62 



- — + ^ + 2-^-^ iky<^ - ikM - -(h'—r—ik 



b a} fc2 



2^fc2 a2"^^^ 



The variation with respect to /32 gives 



2„,' 



2a 



c.c. 



(A5) 



1 



G 



fc262 

k^a"^ 



J 



(A6) 



Here we have defined 



A 



a 
a 



^ 2^ — - 

6 fc,2 a2 6 



After substituting this result into the action, we obtain the action for physical variables: 

^(2)scalar ^ J ^,^^3^. [^GG ^ ^00 ^ ^,7,7 ^ ^G0 ^ ^0,7 ^ ^,7G] 

The term L*^*^ is given by 



(A7) 



(A8) 



--7 |G|^ 
a ' 



4 626'2 fc2 



A2 a2 fc2 \ a 



2a'b 



2u2 



-V 



a^_b^ 
a b 



A2 

bb" + y2 

A 



a 



a b 
a b 



62 pv'^kl fa' 



a^ X k^ \a 6 ' 



1 l/_ 

a b 



\G\' 



Using the canonically normalized variable G = \/2bG, we finally have 



(A9) 



L 



GG 



6 ^ 



2fc2 62 A2 I 



■2a2v^ + :^M_^_^ 
a2 A;2 6 62 



L// L/2 



2a2 /c2 
a b 



4;^^ 1 ( - -^ 



A:2 a2 , 

6^ r 1 

" "6 1 A 



a 
a 6 



|G|2. (AlO) 



The term L'*"*' reads 



2 t'^Jl'^x 



fc2 



2a2 



A2 



2a? fc2 



+ 4 



fc2 
.fc2 



fc2a2 

fc262 



A:2 a2 



Using the canonicaUy normaUzed variable 



A:2a2 ' 
we can write 



^2^,2 ^252 



A kla? 



'|2 



,2 ^ 



A2 



+2- 



2a2 A:2 









fc2 






A;2 


)G 


+ 4^^ ^^'<; 


6' + 


fc,2 a2 





W 2 fc2 52 



a'2 



'fc2 a2 



fc262 
fcp2 



1 

2fe2 



!>20'2 fc252 



Similarly, the term L"^'^ reads 



f2 /r2 

2 fc2l'^ I 



kl\J\' 



2 

/•2l2„,'2 l2 



1,2 J,'2 



k^ 



2a? k^ 

Using the canonicaUy normalized variable J 

k 



Aa^fc| 
2 fc2 



2 



/IfcxI 



J , we obtain 



log 



- log ■ 



4 2 

6 /s A 



A a2 A;,,, 



2a2 /c2 • 



k^ a2 



Aa2 fc| 
2 fc2 



+2 ^ aVy + T^^^/V^ + 4^^ _ ( 2^ - 4^ - 2:^ 



6 



k' X'\ fv' /2 fc^ 



A 



A a2 /ov 



The mixed term L'~^'^ is given by 

&4 



(-2 /2 kx 4>' 

2a^^ kl A 



266' 



,62fc|£ ^a' 
A: ?, A 



/o 2.2T. Pb^v'^kl klb'^b^ 
-frr I 2a^6V + „ , -rf + 4-rf ■ 



A2 
a262 



2a2 /c2 



A:2 a2 



A 



a/_ l/_ 
a b 
62 6' 



66' 



62 /fc2 52^2^,2 ^2 



A2 a2 kl a4 



i-2 



2 62 



" "6 
266' 



fc262 
kla^ 

£62 fc| 
A2^fc2 



{G5(t)* + G* 



Or, using variables G, Scj), we have 



b' 62 fc2 



fc2 A 



6 a2 /c,2 A 



Pv'^ kl 
2a? B 



a" 



a b 
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y (I)' b^ eb^fv'^ki b'(j)'b^ki 



b A2 a2 fc2 



A a2 



A:2 



fc2 



5 A2 a2 



6' 



(G(50* + G*5(j)) 



The terms containing 



and J are given by 



2fc2 



A a2 fc2 J 



626'2 fc2 



J* '(50) 

V' P k^ (j)' k%^ 



2 i-2 



A a2 ky A A:2(j2 



(J(50* + J*6<t)) 



k'^ 

./2 



a2 A:2 



A a2 fc„ 



fv^^Pk, 
2 A a2 A:,, 



(J(50*' + J* (50') 



The above action can be rewritten by using variables 
bL*'' 

kii 



and J 



J as 



2„/2 



0/ 62 fc2 



+ J* (50) 



+ ((.2,2^+^/2.'2§ 



626'2 kl\v'f fc 0' /fc2^,2 



//^z;'2 + a^6V^ 
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~l^'^Yy 
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— — — (J50* 



J 5^* + J* 
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Finally, the terms containing J and G reads 
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kxky 2 I kxky f V 
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Pv''- 
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Using variables G = \/2&G, J 



^%^J, we obtain 



2k^'" 2fc A ' - 
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bb'-^- 



J'G* + J*G 



bb 



A /j,,. 



(JG* + J*G) 



■'^'^ + ir^f'-'^'^T. ) t:t2 f (^G* + rG) + ( ) (^g* + rG) 



p2„y2 "-a: 



+ ( a^b^V + 



2a2- 
2a' 



v' f k 



bb' k fv' f 



2„72 '^x 



J \ a? k. 

k,y. 



y f ky \ X a 

b^b'^ kl\ v^£^l fa' b' 
A ky A 



JG* +J*G) 



(A20) 



Thus, we have obtained the action for the physical variables G, J, and 6(j). 

It is useful to check if J, G and d<f> decouple to each other in the special case ky = for which there exists rotational 
symmetry in y — z plane. Indeed, we find 

^ ^0,7 ^ ^GJ ^ Q (A21) 
The decoupling occurs because the helicity conserves in this special case. The other parts of action become 
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